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A High-Speed Computer for Predicting 
Radioactive Fallout 

J. H. Wright, L. Taback, and H. K. Skramstad 

A very high-speed analog computer has been developed for predicting the geographic 
pattern of radioactive fallout that will occur following a nuclear-weapon explosion. Input 
data consist of wind magnitudes and directions as functions of initial height in the cloud, 
and the distribution of radioactivity in the cloud as a function of the particle size (in terms 
of settling rate) and initial height in the cloud. Output data are presented on a cathode- 
ray tube, so that the luminance at any point on the screen represents the intensity of fallout 
accumulated at that point. 



1. Introduction 

The computer described in this paper was de- 
veloped for the U. S. Weather Bureau and the 
Albuquerque Operations Office of the Atomic 
Energy Commission, with the immediate objective 
of operation in the Pacific Proving Grounds during 
the 1956 tests of nuclear weapons. Two essentially 
identical units were built and used by the fallout- 
predict ion group as aids for safety purposes. These 
units, as reported here, were experimental models 
lacking engineering refinements. A quite different 
type of special-purpose analog computer, developed 
by the Sandia Corp., was also employed. Another 
machine, using mechanical-optical techniques, is 
under development at the Los Alamos Scientific 
Laboratory. The most distinctive features of the 
National Bureau of Standards machine are its 
virtually instantaneous speed of solution, straight- 
forward procedure for setting in data, and adapt- 
ability to variations of the "model" of the radio- 
active cloud. Various limitations, as well as some 
possible refinements, are discussed. 

In the past, hand methods of solution have been 
used extensively. These require considerable sim- 
plification of the problem and exploitation of the 
operator's judgment to achieve the solution in a 
reasonable time, e. g., an hour or less. Very large 
digital computing facilities have also been used 
with considerable reduction in the computation 
time. Both of these techniques will be of continuing 
importance; the hand methods, because they re- 
quire negligible outlay for equipment and installa- 
tion; the large digital computers, because they can 
do increasingly complicated problems simply by 
taking more time. In contrast, the desired objec- 
tives for the special-purpose computer were: 

1. Small size and reasonable portability. 

2. Low cost— under $10,000. 

3. Problems set up without the need for mathe- 
matical training in numerical analysis or digital 
programing. Actually, the operator need not be 
familiar even with mathematical functions such as 
cosine or exponential. 

4. Very rapid solution of the problem. 

5. Direct visual presentation of the final solution. 

6. Rapid problem modification. 



These requirements have generally been met in 
the first model, with some reservations. The solu- 
tion time is considered to be "instantaneous" be- 
cause it is negligible compared to the time required 
for human observation of the output display of 
data. The estimate of cost is crude and is based on 
construction of 3 or 4 duplicate machines, exclusive 
of engineering time or test or of commercial profit. 
The over-all accuracy of the machine has not been 
studied systematically, because of the immediate 
need for field service, and cannot be discussed. The 
most serious source of error was the conventional 
cathode-ray tube employed, which, however, pro- 
vided the most appealing and useful feature of its 
operation in the first application. 

It should be noted at the outset that this computer 
is not a direct simulation type of analog machine, 
but rather a scheme for using analog techniques to 
mechanize the problem in a maimer somewhat 
similar to the organization of the problem in a large 
digital computer. In terms of apparatus, this 
machine consists primarily of function generators 
and scanning voltages. Time, as represented in the 
flow of machine operations, has no relation to the 
time variable in the original physical model, but is 
used for sequencing purposes only. 

2. The Fallout Problem 

The assumptions made in the solution of the 
problem are not unusual in the field of fallout cal- 
culations, and will be listed here in some detail. 
First, it is assumed that a cloud of radioactive 
particles has formed, with a definite shape and size, 
and, at some instant considered as zero time for the 
calculations, the various kinds of particles at various 
locations in the cloud begin to settle, being displaced 
laterally by the winds as they fall to the ground. 
This ignores the possibility that some particles, 
notably very large ones, may have reached their 
maximum height and started falling while other 
types of dust particles are still rising toward their 
"initial" positions. The cloud is assumed to have a 
vertical axis of symmetry, whereas actually the 
horizontal winds can tilt and otherwise distort the 
cloud while it is forming. Vertical winds are neg- 
lected. The large-scale vortex motions caused by 
the forces developing the cloud, and any effects on 
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the wind pattern resulting from the energy of the 
burst, are also neglected. It is assumed that the 
characteristics of the particles are sufficiently 
described by specifying a rate of fall as a function of 
altitude, and a particular weighted value of radio- 
activity for each of the various ranges of fall rate 
and initial height. We will hereafter refer to a 
"slowness" variable, ,9, because the corresponding 
variable in the machine is proportional to fall time 
rather than rate. The slowness varies greatly with 
particle size. However, it is assumed that the 
relative times spent in the various altitude layers are 
the same for all particles starting from the same 
height. Thus, the law of fall can be made correct 
for some intermediate size of particle, but becomes 
increasingly incorrect toward the extremes of particle 
size. 

Some of the above assumptions become increas- 
ingly important for calculations at positions nearing 
ground zero, the site of the explosion. Eventually, 
perhaps within an area of two mushroom diameters 
centered on ground zero, these factors pretty well 
invalidate the detailed computations. The mush- 
room diameter refers to the mean diameter of the 
widest portion of the particular cloud considered. 

A basic assumption is that the wind at each level 
is invariant over the entire area of fallout, and does 
not change during the entire time of fallout. Time 
and space variability of winds therefore introduce 
increasingly large errors for particles taking a long 
time to fall, and for fallout far from ground zero. 
The relative importance of such variations during 
actual operations is the subject of present research 
by Sherman, 1 Nagler, 2 and others. Detailed calcu- 
lations of the effects of variability at all altitudes 
and areas, and during the time of fallout, could 
complicate the problem enormously — by several 
orders of magnitude. Thus, practical considerations 
of operating conditions become all important, as to 
whether the requisite detailed weather data can 
actually be obtained, and soon enough to be valid. 
For example, if a fallout prediction is to be made 12 
hr before detonation of a bomb, and must allow for 
particles requiring 24 hr to reach the ground, a very 
detailed 36-hr weather forecast is necessary. Such 
forecasting is difficult, if not impossible. Roughly 
speaking, we may say that beyond about a 200-mile 
radius, or 12-hr falltime, fallout predictions neglect- 
ing variability are of doubtful validity, at least 
without extensive reappraisal by weather experts. 

The radial distribution of radioactivity at any one 
initial height in the cloud is assumed to be the same 
for all types of particles. As will be seen, the present 
computer does not control this radial distribution 
closely. Thus, it may be roughly uniform, or fall 
off with increasing radius, depending on whether or 
not a compensating waveform is applied to the in- 
tensity grid of the cathode-ray tube. (An arbitrary 
profile could be generated, but such a refinement 
was considered unimportant in the first model.) 
Finally, it is assumed that no diffusion occurs during 
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fallout, so that particles of a given size originating 
in any specified layer of the cloud maintain their 
relative positions as they fall. Thus, the calculation 
of wind effects on particles on the axis is sufficient; 
i. e., the entire "slice" of cloud (considering particles 
of one size range) moves as a piece during the fallout. 

The machine does not need to distinguish between 
the number of particles per unit volume, p i} and the 
activity per particle, q iy for particles of a given size 
(or "slowness"), i. The machine input known as 
activity, a, is the weighted activity for that size of 
particle for an entire layer of cloud. Thus, a tj = 
qtijPijVj, where the subscript j refers to the jth layer 
of cloud. Vj is the volume of this layer, so that 
PijVj is the total number of particles in the jth. layer 
of cloud. For a given type of nuclear weapon, values 
must be assigned to the various q ih p tj andT 7 ^, so as to 
define a function of two variables, a=f(h,s), where h 
represents height in the initial cloud, and s the slow- 
ness (time to fall a unit height) typical of a given class 
of particle. This function a is commonly called the 
"cloud activity model." In practice, the cloud 
activity model can be supplied to the machine 
operator as a family of curves, and he need not even 
be familiar with mathematical functions. The 
activity function is set up so that the output is in 
terms of percentage of total activity in the cloud. 
The family of curves are adjusted with h as a parame- 
ter for 0, 20, 40, 60, 80, and 100 percent 3 of the 
maximum altitude of the cloud. This normalizing 
procedure provides for solutions over a considerable 
range of weapon sizes without readjustment of the 
cloud activity model as long as the shapes of the 
clouds remain the same. 

It must be noted that the machine does not make 
any detailed allowance for decay in radioactivity, 
either during the time of fall or afterward. The 
assignment of a tj to a given range of particles at a 
given initial height must include the allowance for 
decay. Furthermore, for a particular activity model 
set up in the machine, the operator cannot make a 
simple single adjustment to show the fallout at some 
arbitrary earlier stage, i. e., before all the particles 
included in the model have reached the ground. 
However, the activity model can be modified to 
include only those particles that will reach the ground 
before the specified time. In any case, the model 
can be set up on either of the following bases: 

Integrated dose. The total dosage from time of 
deposition to some specified time, such as 48 hr., or 
10 years, after the burst. Often the "infinite dose," 
for complete decay, is used in estimating possibility 
of hazard to permanent residents in an area. The 
unit of measure is the roentgen, R, or milliroentgen, 
mr. 

Dose rate. The intensity of radiation at some 
selected time after fallout, such as 48 hr., This 
quantity may be expressed in roentgens or milli- 
roentgens per hour, and is appropriate to considera- 
tion of evacuating personnel or exploring a con- 
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total time of fall, and will not be exactly 20, 40, 60, and 80 percent of the cloud 
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taminated area. It is a little more 4 convenient than 
the integrated dose in deriving a model, because the 
time of arrival need not be known. In practice, it is 
common to use one law of decay applying to all the 
batches of particles, the gross fission-product decay 
expression, A=A t~ 1 - 2 , in which A is the dose rate 
measured at time t after the explosion. A Q is the 
dose rate at unit time, such as roentgens per hour at 

1 hr. On this basis, the completed fallout calcula- 
tions for dose rates, at some time such as 48 hr. after 
the explosion, can be converted by a single factor to 
provide the dose rates at some other time, such as 
10 days later. 

Despite this simplification of the problem, the 
required amount of computation for a nearly in- 
stantaneous presentation is quite formidable. To 
illustrate this, let us consider the following break- 
down of the problem into discrete elements appro- 
priate to digital computation. The cloud is divided 
into 20 layers vertically, and 25 elementary cross- 
sectional areas in each layer. Thus, there are 500 
elemental volumes in the cloud, each containing a 
distribution of particle sizes, with attendant varia- 
tions in "slowness" of fall. If we divide this range 
of particles into 10 classes, there are effectively 5,000 
lumps of material in the cloud, each requiring calcu- 
lation of displacement by the wind in each wind 
layer through which it falls. On the average, 10 
multiplications must be performed to obtain the final 
horizontal displacement of each elementary lump of 
material as it falls to the ground, each the product 
of a time of fall and a wind vector. If the solution 
is to be completed within 1/10 sec, the required com- 
puting rates are: 

(a) 50,000 function table values per second, for 
radioactivity. 

(b) 500,000 simple multiplications per second, for 
partial displac erne i i t s . 

(c) Appropriate combination of these partial 
results, and presentation on a map. 

The display device would be required to generate 
"picture elements" of varying intensity at the rate 
of 50,000/sec, one for each representative lump of 
radioactive dust at its final fallout location. These 
picture elements tend to be superimposed, because 
light particles originating in a given layer will fall 
on top of heavier particles originating from lower 
altitudes in positions relatively "downwind" with 
respect to the light particles. This situation will be 
referred to as overlap. Thus, the display device 
must be able to accumlate these overlapping contribu- 
tions, if the computer proper is to be relieved of this 
burden. 

The accuracy requirements on individual opera- 
tions are rather low, generally ranging from errors of 

2 to 10 percent. An analog technique that endeavors 
to meet approximately the requirements of the 
example will now be described. 

3. Wind Computations 

Let us assume that the mushroom cloud (see fig. 1) 
is divided into horizontal layers at graded intervals of 




-r*-kh = I 
Figure 1. Mushroom cloud. 




Figure 2. Wind hodograph. 



altitude such that a particle of typical size spends 
equal times falling through successive layers. Let 
h represent the initial height of any particle. Now 
let us define a variable A*, a function of h, which is 
proportional to the time for a particle to fall from 
any specified height to the ground. The variable 
A* has the value 1 at the top of the bottom layer, and 
n at the top of the nth layer. Let V t represent the 
average speed and 6 t the average wind direction in 
the ith layer, 6 being measured clockwise from north. 
The cloud diameter as a function of height is repre- 
sented by d(h). 

Figure 2 shows a wind hodograph. The dotted 
line may be considered to represent the horizontal 
projection of the path of a hypothetical weather 
balloon having a constant rate of ascent. If an 
appropriate altitude scale is marked along such a 
curve, a vector such as T 2 joining altitude points 
hi and h 2 would represent the direction and magnitude 
of the average horizontal wind exerted on this 
hypothetical balloon as it rises completely through 
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the layer between heights hi and h 2 . Such hodo- 
graphs are commonly used to obtain measurements 
of winds aloft. In the case of the descending motion 
of particles, the hodograph would be drawn in the 
reverse sequence, starting with V 6 instead of \ x . 
However, the total relative displacement would be 
the same, if the time spent in each layer were the 
same as for the balloon. The horizontal displace- 
ment of a particle while it falls from h 2 to hi is 
obtained by multiplying this wind vector by the 
time required for this particular particle to fall 
through this layer. A particle of size a falling from 
the top of the third layer to the ground would be 
displaced in an east-west direction by the amount 

x=s(a)[Vi sin 6±-\-V 2 sin 2 +T^3 sin 3 ], 

where s(a) is the slowness for this size of particle, 
the time to fall any one height interval. A similar 
cosine expression holds for the north-south compo- 
nent of displacement. 

It is apparent that there is an approximation 
introduced by assuming a finite number of wind 
layers, because the dotted line of the wind hodograph 
deviates from the straight-line vectors joining the 
layer boundar}^ points h u h 2 , etc. However, this 
introduces no error for the particles falling com- 
pletely through a given layer, because the total effect 
of the wind is specified correctly. There is a small 
error, within the layer only, for those particles that 
start to fall from some height within the particular 
layer under consideration, because the analog com- 
puter interpolates along the straight line in space. 
If a digital computer were used, it would be necessary 
to divide the cloud into a great number of layers, 
or to provide an interpolation routine. 

The basic equations for the final ground position 
of a particle of slowness s j} originating on the axis 
of the cloud within a layer for which n<^h*<^n-\-lj 
are therefore 

x=s j ibVt sin 6i+(h*—n)V n +i sin 6 n+l 



y~- 



r n 
=*i J2V f cos di + (h*-n)V n+1 cos 6> 

Li=l 



+1 \> 



where x represents the fall position east of the 
explosion site and y the position north of the explo- 
sion site. In each equation, the second term., such 
as (h*—n)V n+ i sin0 w+1 , represents the initial displace- 
ment of the particle while falling through a fraction 
of layer (n+1) enroute to the top of laj^er n. 

4. Description of the Computer 

A simple block diagram of the apparatus is shown 
in figure 3. From the viewpoint of the operator, 
there are four basic units: 

(a) The winds unit, with knobs for setting in the 
magnitudes and directions of the mean winds in 
each of the 20 layers. 
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Figuke 3. Block diagram of computer. 



(b) The cloud geometry unit, including the 
cloud diameter and cloud raster functions. There 
are several knobs for adjusting the cloud shape, and 
one knob for setting in an over-all scale factor for 
cloud size. 

(c) The cloud activity unit, which has 36 knobs 
for setting in the "model" of distribution of radio- 
activity in the cloud. 

(d) The display unit, which has x and y scale-factor 
controls (marked y-eanp and x-amp), a brightness- 
level control, and a brightness-slope control. 

There is also a central control unit containing all 
timing and sweep circuits, summing amplifiers, 
slowness modulators, and final x and y driver 
amplifiers. It needs no attention in normal opera- 
tion. In addition, there is a photometer used to 
measure the brightness on the screen as a direct 
measure of fallout. It is an electrically and mechan- 
ically separate entity. 

For explanation, let us consider first that the 
cloud is collapsed laterally so that all the particles 
are on the vertical axis of the cloud. Furthermore, 
consider that we are dealing only with the lightest 
(i. e., slowest) particles given any consideration in the 
problem. We will later modify the solution to include 
heavier particles, and those off the axis of the cloud. 
A periodic sawtooth voltage is generated to represent 
the variable A*. This is called the height sweep. 
The voltage rises linearly with time to a nominal 
+ 100 v, is reset to zero in about 40 ptsec, and repeats 
each 1/10 sec, scanning through the full range of A* 
each time. This is the basic cycle of computer 
operation, and the entire solution is generated anew 
with each sweep through the full cloud height. This 
voltage is applied to the winds unit, a function 
generator with 20 elements sequentially activated 
as the input voltage rises, but with their outputs 
added together at all times. At the start, with the 
A* voltage at 0, none of the elements produces an 
output. As the A* voltage rises from to a nominal 
value of +5 v, the output of the No. 1 element, 
representing the lowest wind layer, also rises from 
to +5 v. As A* rises from +5 to +10 v, the 
output of element 2 becomes activated, and rises 
from to +5 v, thereafter remaining at +5. This 
goes on until A* reaches a nominal value of +100 v, 
and all 20 function-generator elements are pro- 
ducing outputs. 
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Figure 4. Schematic of winds unit. 

This sequential and cumulative action provides 
the necessary effect of only 1 wind contributing to 
the total displacement of particles falling from the 
top of the first layer, but 20 winds contributing 
to the total effect on particles originating at the 
very top of the cloud. For a particle originating, 
say, halfway in layer n, the cumulative wind effect 
is represented by the sum of the earlier sustained 
voltages, representing lower altitude winds, plus 
half the full scale voltage from the partly activated 
nth element. The situation is shown in figure 4 for 
particles in the first five layers of altitude. Actually, 
before these voltages are summed, they pass through 
attenuating potentiometers, adjusted according to 
the mean magnitude of the wind in each layer. The 
nominal +5 v represents a maximum-scale wind, 
e. g., 100 knots. If there is no wind in a given layer, 
that potentiometer is simply set to zero. The 
voltages are also resolved into x and y components 
by sine-cosine potentiometers set according to the 
wind direction, one for each layer. 

Figure 4 also provides a simplified electrical dia- 
gram of one element, for one wind layer. When A* 
is at the beginning of the sweep, the lower diode, 
connected to a nominal +5 v (about 4 v is used, 
to allow for forward drop in the diodes), is non- 
conducting because it is held in a reverse bias con- 
dition by the action of the large negative voltage 
— V. The upper diode is conducting, meanwhile, 
and sets the input to the wind potentiometer at v. 
The anode of this diode is shown connected to +0.6 
v, a value adjusted in operation to compensate for 
the diode forward drop so that the wind-potentiom- 
eter output actually is 0. As the A* voltage rises, 
the voltage at the arm of the altitude breakpoint 
potentiometer will at some time, depending on its 
setting, rise above v. While it is between and 
+ 5 v, both diodes are nonconducting, and a voltage 
builds up across the wind potentiometer, reaching a 
nominal value of +5 v before being limited by con- 
duction in the lower diode. 

The output of the wind potentiometer drives one 
side of a sine-cosine potentiometer. An intermediate 



buffer amplifier is necessary to prevent electrical 
loading of the wind potentiometer, and a 1:1 inverter 
amplifier is necessary to drive the other side of the 
sine-cosine potentiometer. i?i is necessary for elec- 
trical reasons to prevent excessive current being 
drawn from the breakpoint potentiometer at high 
settings, and to standardize the forward current 
through the diodes. Resistors 7? 2 are identical 
summing resistors feeding ordinary analog summing 
amplifiers, a total of 20 connected to the east-west 
summer, and 20 to the north-south summer. The 
altitude breakpoints are set to occur at 20 steps of 
the A* voltage, equally separated from to +95 v 
(for an A* amplitude of 100 v). The value of the 
nominal +5 v supply is adjusted to make the indi- 
vidual steps add up properly. The 20 rheostats 
labeled i? 3 are calibration adjustments to correct for 
component variations, thus making the full-scale 
contributions from all the wind potentiometers equal. 

As front panel adjustments, the R% rheostats could 
be used to weight the winds so that the 20 layers are 
equally spaced heights, rather than irregular height 
intervals spaced at equal intervals of fall time (for a 
selected particle size). This may prove to be de- 
sirable for convenience and flexibility. If this is 
done, the settings must be changed when shifting to 
a problem with quite different initial cloud height, or 
choice of typical particle. 

Additional sawtooth voltages are generated by the 
particle slowness modulators. Each modulator pro- 
duces four sawteeth while the height sweep scans 
through each cloud layer. The amplitudes of these 
sweeps are scaled according to the relative time of 
fall. At the midpoint of the sweep, for example, the 
vector sum of the outputs of the x and y channels is 
half the full amount. The result is to trace a radial 
line, or "spoke," that is the locus of fallout of 
particles of various kinds up to the slowest-falling 
type. These spokes are nominally constant-height 
lines, but are curved somewhat, because there is a 
little variation in height, and hence in mean wind, 
during the scan of particle slowness. The effect is to 
multiply each component by the sawtooth voltage, 
but no conventional multipliers are used. The effect 
is achieved by use of voltage-controlled sweep ampli- 
tudes, of whatever polarity is dictated by the wind 
pattern. These modulators are actually integrating* 
circuits, using operational amplifiers. They are 
reset 800 times per second by triode reset circuits 
that require about 20 /xsec. 

To take into account the diameter of the cloud, the 
spot on the cathode-ray tube must be effectively 
enlarged to have a diameter proportional to that of 
the cloud, changing with the voltage A*, as the solu- 
tion of the problem progresses. To spread the spot 
sufficiently, it has been necessary to resort to a high- 
frequency spiral raster, with a circumferential fre- 
quency of about 100 kc and a repetition rate (raster 
rate) that varies from 4 to 20 kc or more. The high- 
est frequency is associated with the smallest spiral 
diameter, because the spiral advances outward at a 
fixed rate until it reaches a specified voltage repre- 
senting the cloud diameter. Thus, larger raster 
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diameters will be repeated at a lesser rate, and will 
have more turns in the spiral. The computer unit is 
called the cloud raster generator. It is driven by a 
voltage from the cloud diameter generator, which is 
a conventional diode function generator approxi- 
mating a function of one variable (the A* voltage) by 
means of connected line segments having adjustable 
slopes and breakpoints. The x and y voltage com- 
ponents of the spiral are inserted into the final x and 
y driver amplifiers as increments to the basic fallout 
positioning voltages. 

The cloud activity generator receives the A* and s 
voltages and continuously generates an output that 
is a function of these two variables. The output of 
this unit can be displayed on the main cathode-ray 
tube when desired, in the form of a family of curves 
of activity a versus slowness with A* as a parameter. 
Actually, 80 curves appear in the family, but the 
operator sets in curves only at the top of the cloud, 
ground level, and four intermediate points at equal 
intervals of A*. Each curve is made up of five line 
segments. This generator must be unusually fast in 
operation, because s sweeps through its entire range 
of 50 v in 1/800 sec. Furthermore, it should be possi- 
ble to adjust to arbitrary functional variations be- 
tween curves in the family, rather than use poly- 
nomial methods. It is desirable that the operator 
simply adjust a row of knobs until the proper func- 
tion of activity versus slowness is obtained at ground 
level, then independently adjust to match a specified 
curve at the next selected height in the cloud, and so 
on to the top of the cloud. 

Function generators developed by the George A. 
Philbrick Co. were emplo}^ed in the first two com- 
puters. They permit such independent arbitrary 
adjustment. However, the interpolation between 
curves for a constant value of the s variable is not 
linear. For example, in progressing from a hori- 
zontal line to one with a V-shaped notch, the inter- 
mediate curves are truncated, with flat sections, 
rather than V -shapes of intermediate height. The 
device is all electronic, using a matrix of selector 
diodes. Better interpolation is desired, and an 
electronic function generator 4 recently developed by 
L. Taback of the Bureau may find application in later 
models of this computer. The Taback generator 
uses a time-division technique for interpolation be- 
tween curves in a family. 

The output voltage of the activity generator is 
connected to the intensity grid of the final display 
cathode-ray tube. It thereby modulates the elec- 
tron beam so as to make the brightness proportional 
to the activity of the particle type at the initial height 
under consideration at any instant. The intensity 
level of the cathode-ray tube does not remain 
sufficiently constant (at a specified grid test voltage) 
for a permanent type of calibration. Furthermore, 
the luminance versus grid voltage curve departs 
considerably from linearity in some regions of opera- 
tion. Hence, slope and brightness-level controls 
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Figure 5. Basic waveforms 
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* L. Taback, A function generator for two independent variables, Control Eng. 
3, 129 (1956). 



Figure 6. Output of north-south winds summer. 

are provided. The slope control is an attenuator 
in the activity signal line. These controls must be 
adjusted with the entire computer in operation after 
insertion of a calibration problem with known winds, 
cloud geometry, and activity model. The procedure 
is to first set the x and y deflection factors for the 
cathode-ray display so as to make the display read 
correctly in miles for the eastward and northward 
limits of the pattern. For this purpose the test 
problem can be made very simple. Then the over-all 
size factor for the cloud is adjusted to make the 
mushroom diameter appear to be the proper num- 
ber of miles on the screen. The slope and bright- 
ness-level controls for the cathode-ray tube are then 
adjusted to make the photometer read as correctly 
as possible in areas with various ranges of brightness. 

Some of the basic waveforms in the computer are 
shown in figure 5. All are exhibited on the same 
time scale synchronous with the A* voltage, but to 
simplify the diagram, only 5 wind channels are shown 
instead of the actual 20, and only 10 s sweeps instead 
of the actual 80 during each A* sweep. Four of the 
waveforms are associated only with wind 2, while 
the two waveforms at the bottom show the contribu- 
tions from channels 1,3,4, and 5 as well. 

Figures 6 to 14 were taken before the first com- 
puter was completed. They do not exactly repre- 
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Figure 7. Output of north-south slowness modulator. 
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Figure 8. Output of activity generator 



sent normal operation, providing less than the nor- 
mal 80 "spokes," for example. Because it was not 
possible to take pictures subsequently, these are 
presented to aid the discussion. 

Figure 6 is an oscillograph trace 4 of the output of 
the north-south summer versus A* voltage. Figure 
7 then shows the result after the waveform shown 
in figure 6 is used to modulate the particle-slowness 
voltage. The triangular waveforms actually have an 
envelope that follows the waveform of figure 6. 

Figure 8 is an oscillograph picture of the output 
of the activity generator, set for convenience of 
display, and not to represent any actual explosion 
cloud. The abscissas represent particle slowness, 
and the ordinates represent the output activity, a. 
Each of the individual curves corresponds to a differ- 
ent value of A*. 

Figure 9 shows the display of fallout on the cath- 
ode-ray-tube screen in the absence of any intensity 
modulation by the activity generator, and without 
the spreading of the beam by the cloud raster genera- 
tor. Each spoke is the result of one slowness sweep. 
Fallout from successively higher initial layers in 
the cloud corresponds to progression of the spokes 
in a clockwise sense for this case. The number of 
spokes is determined by the ratio of the s sweep 
frequency to the A* sweep frequency. 

Figure 10 was photographed under the same com- 
puter settings but includes the intensity modula- 
tion according to the function a. 
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Figure 9. Simple spoke display. 




Figure 10. Intensity-modulated spokes. 

Figure 11 is similar to figure 10, except for more 
drastic modulation of the intensity. Here only the 
tips of the spokes show, as would be the case if 
only the slowest particles were radioactive and 
numerous enough to be of importance. 

Figure 12 shows the result of superposing the 
cloud raster on the display of figure 11. In actual 
use, the oscilloscope spot is defocused so that the 
spiral pattern does not show clearly. (The appear- 
ance 4 of spirals of widely different diameters on the 
same centers is not proper, and is the result of 
faulty synchronization existing at the time of this 
early picture.) 
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Figure 11. Spoke tips 
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Figure 12. Spiral rasters. 

Figure 13 shows how the final display appears with 
all units of the computer in action. The computer 
was adjusted to produce this apparent fallout with- 
out knowledge of the characteristics of any actual 
bomb, hence the pattern is not representative. 

Figure 14 shows one of the computers. The 
cabinet on the right contains the computer proper, 
together with an auxiliary 5-in. oscilloscope. The 4 
chassis below the 5-in. scope form the winds unit, 
each containing 5 wind channels. Below these, mov- 
ing downward, are the cloud activity unit, central 
control unit, cloud diameter unit, and cloud raster 




Figure 13. Complete fallout display. 




Figure 14. Computer. 

generator. The 21 -in. cathode-ray display unit is 
mounted on a table at the right. On top of it is the 
meter of the photometer, a modified commercial unit 
with a probe area % in. in diameter. The probe is at 
the end of a cable. The cabinet at the left contains 
power supplies, test equipment, and a spare chassis. 
The second large cabinet is unnecessary, and the 
essential power supplies can be mounted in a small 
unit forming a pedestal for the display unit. Forced- 
air cooling is employed, with lateral airflow from 
front to rear over the chassis. The total a-c power 
requirement is about l,500w. 



108 



5. Comments on Operation 

Because the machine is used intermittently, and 
accuracy requirements are relatively low, chopper 
stabilization was not employed in the operational 
amplifiers. Thus, balance adjustments are neces- 
sary. These are made after 1 hr. warmup time, the 
process taking about 5 min. Following this, some 
readjustment is necessary at intervals of 2 to 6 hr. 
Checks on the brightness level of the scope and 
balance of the activity generator are necessary at 
about the same intervals. A few other basic elec- 
trical adjustments are made at intervals of a few 
days, including sweep amplitudes and altitude 
breakpoints. 

In use, provided the cloud geometry and activity 
model had already been inserted in the machine, a 
completely new pattern of winds could be inserted in 
about 2 min. once winds were specified at the proper 
altitudes. Results of problem modification were im- 
mediately observable on the display, but 2 to 5 min. 
were required to take a useful pattern of measure- 
ments with the photometer probe. 

Modification of the activity model from a set of 
tabular data could be done within 5 min. The shape 
of the cloud could be altered as desired, but it seldom 
needs serious adjustment. The over-all size of the 
cloud could be altered instantly by turning a single 
knob. 

Relative position of ground zero on the screen 
could be adjusted as quickly as desired, but with 
some limitations due to shortcomings of the oscillo- 
scope. Likewise, the picture could not be expanded 
greatly, because of the limitation in cathode-ray 
deflection speed to about 0.6 in./^sec. This meant 
that the mushroom could not appear as more than 
1% in. in diameter on the screen, due to the periph- 
eral rate of 100,000 turns/sec in the spiral. An 
overlay map on transparent backing was used over 
the display screen. 

6. Cathode-Ray Tube Problems 

Because of "overlap," any elementary area on the 
cathode-ray screen may receive many contributions 
of beam current during any one computation cycle. 
After one "strike" by the beam, others will occur as 
the spiral raster moves on out along a given spoke, 
or along an adjacent spoke. It might seem desirable 
to use a phosphor screen having long persistence, so 
that it would integrate these contributions. How- 
ever, the accumulated light-energy output versus 
electrical charge resulting from the various strikes is 
not a correct integral. The relation between beam 
current, exposure time, and previous history is quite 
complicated and unsatisfactory. Storage tubes and 



dark trace tubes were 4 considered, but the practical 
choice was a short-persistence 4 tube. Ideally, the 
persistence should be negligible. The summation is 
accomplished in the averaging action of the pho- 
tometer deflection. Because visual display was 
desired, the computer repetition rate had to be fast 
enough to reduce the flicker in visual effect. A 
compromise of 10 cps was made in this computer, 
although 30 cps would be necessary for really sub- 
stantial reduction of flicker. Actually, the phosphor 
was a fairly severe compromise, being a P-4 type. 
Theoretically, a considerably faster type would be 
desirable. 

7. Refinements 

Fairly obvious engineering refinements will reduce 
size and improve stability, particularly in choice of 
components. The number of components can be 
reduced, and the winds unit, in particular, could be 
transistorized. Zener voltage diodes can be used to 
simplify the winds unit and eliminate the low- 
voltage precision power supply. 

Greater accuracy and usefulness under some con- 
ditions could be achieved by eliminating the cathode- 
ray tube altogether. The same basic circuits would 
be used, but could be run much slower, thus reducing 
the basic problem of achieving a sufficiently high 
raster rate relative to the slowness sweep. The out- 
put of the activity generator would be gated into an 
averaging circuit when and only when x and y simul- 
taneously have values falling within a specified 
elementary geographical area. For any such elec- 
trically probed area, the gate would open many 
times during each solution time, due to the overlap- 
ping previously discussed. The output would be 
read directly on a meter, or could be printed out on a 
standard electric typewriter operating from a rela- 
tively slow analog-digital converter. The presenta- 
tion could have the form of a map, each advance of 
the typewriter representing a standard number of 
miles in a routine scan of the fallout area. 

The present computer design could be extended to 
accommodate more elaborate formulation of the 
problem, such as time or space variability of winds. 
However, it is not known whether it is entirely de- 
sirable to attempt this. 



The authors are indebted to Lester Machta of the 
U. S. Weather Bureau, who initiated the project on 
behalf of the Atomic Energy Commission, and to 
K. M. Nagler, also of the U. S. Weather Bureau, for 
his assistance during the project. 

Washington, September 13, 1956. 
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